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We have applied the variational ii-matrix method to calculate the reflection and tunneling prob- 
abilities of particles tunneling through one-dimensional potential barriers for five different types 
of potential profiles - truncated linear step, truncated exponential step, truncated parabolic, bell- 
shaped, and Eckart. Our variational results for the transmission and reflection coefficients are 
compared with exact analytical results and results obtained from other numerical methods. We 
find that our results are in good agreement with them. We conclude that the variational ii-matrix 
method is a simple, non-iterative, and effective method to solve one-dimensional quantum tunneling 
problems. 
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I. INTRODUCTION 

Quantum tunneling is of considerable interest in mod- 
ern physics [1]. Tunneling of particles through poten- 
tial barriers have applications in different areas such as 
radioactive decay [2], black hole production [3], optical 
wave-guides, lasers, quantum well structures [4], and elec- 
tron transport in micro and nanoscale devices [5] . In fact 
spin polarized tunneling in semiconductor and magnetic 
heterostructures has opened up a new area in condensed 
matter physics, known as "spintronics" [6]. 

The tunneling coefficient is the key element of the tun- 
neling process that must be understood. For example, in 
the case of nano-devices it determines the current-voltage 
characteristics. In a tunneling problem one usually solves 
for the tunneling coefficient by solving the Schrodinger 
equation with appropriate boundary conditions. But 
there are only a few potentials for which the Schrodinger 
equation can be solved analytically so that exact results 
for the tunneling coefficients can be obtained. Other po- 
tentials must be solved numerically using techniques such 
as the Wentzel-Kramos-Brillouin (WKB) [2], modified 
conventional WKB (MWKB) [7], modified Airy function 
(MAF) [8], numerical matrix (MAT) [9], transfer matrix 
(TM) [10], and E-matrix [11] methods. All these tech- 
niques have been employed to calculate the transmission 
coefficients for different types of potential barriers. 

The i?-matrix method was originally developed by 
Wigner for use in nuclear physics [12]. Later his method 
was extended by others to solve a variety of different 
problems in atomic and molecular physics [13-16]. A 
variant of the 7?-matrix method known as the variational 
i?-matrix method, has been successfully applied recently 
to the problem of quantum tunneling through a potential 
barrier [17]. 

In this paper, we report on our test of the effective- 
ness of the variational i?-matrix method in obtaining the 



transmission coefficients for a variety of one-dimensional 
potentials of interest. We find our variational results to 
be in good agreement with existing analytical and nu- 
merical results. 

In Sec. II we briefly present the variational i?-matrix 
method. In Sec. Ill we apply the method for particle tun- 
neling through potential barriers. Our results are com- 
pared to results from analytical expressions and other nu- 
merical methods. The results are summarized in Sec. IV. 



II. FORMALISM 

In this section, we briefly outline the variational R- 
matrix method following Ref. [17]. We begin with the 
one-dimensional time independent Schrodinger equation 
for a particle of mass m and with energy E in a potential 
V(x) (assumed to be real), 

+ - (£ _ ni)1 ^„. (1) 

As Eq. 1 is a linear and second order differential equation 
there are two independent solutions, tp and "0 with the 
same energy E. We define three regions: an inner region 
(a < x < b) where the potential is a function of the 
position x, and two outer regions defined by x < a and 
x > b, where a and b are the boundaries separating the 
inner region II from the outer regions I and III. 

In the regions I and III the potentials V\ and V 3 are 
either zero or constant. So the wave-functions (?/>i, f/'i) 
and (ip3, ips) are oscillatory and can be written in terms 
of exponential functions, 

ipi{x) = a\ exp(ikix) + b\ exp(— ik%x), (2) 

ipi{x) = 5,1 exp(ifcia;) + b\ exp(— ik%x), (3) 

tpa(x) = a 3 exp(ik 3 x) + b 3 exp(-ik 3 x), (4) 

ip 3 (x) = a 3 exp(ik 3 x) + b 3 exp(-ik 3 x), (5) 
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where fa = v / 2m(E - Wj/K and k 3 = ^2m{E - V 3 )/Ti. 
We assume E > V\, V 3 , so fa and £3 are real. The coeffi- 
cients ai,ai,bi,bi (i = 1,3) can be calculated by match- 
ing the wave-functions and their first derivatives at the 
boundaries, i.e. at x — a and x — b. 

In region II, we expand the wave-function in terms of 
a suitable set of N basis functions Xi( x ), 



N 

i>2{x) = ^2c lXl (x). 

i=l 



(6) 



Here c,'s are the variational coefficients. In all the cases 
we have studied we have used Xi( x ) — cos(nix) to be a 
simple basis set. The parameters Kj's are chosen by trial 
and error so that the basis functions \% ( x ) are sufficiently 
complete for our input energy values. 

We set-up the i?-matrix for region II by multiplying 
both sides of Eq. 1 with ij)2{x) and then integrate by 
parts to obtain, 

- / taW 2 (x)] 2 + -p / dx[E-V(x)]$(x) 

+ X(b)^ 2 2 (b) - X(a)^ 2 (a) = 0. (7) 

In Eq. 7 we have defined the logarithmic derivatives of 
■02 at x = a, b: A(a) = tp' 2 (a)/ip 2 (a), A(6) = %}) 2 {b) / ip 2 (b) . 
We substitute Eq. 6 in Eq. 7, which now takes the form 

N 

<3= E [AijW + KV&ijCiCj - A(a)A«. CiCj ] = 0, (8) 

where A, A a , and A b are NxN real symmetric matrices 
defined as 



dx x'i(x)x'j(x) + 



2mE 



dx Xi(x)Xj{x) 



2m 



[ dx Xi(x)V(x)xj{x) 

J a 



= Xi(a)Xj(a), 
^ = Xi(b)Xi(P). 



(9) 
(10) 
(11) 



Q in Eq. 8 does not change by varying a's. Thus by 
taking the derivative of Q with respect to Ci and setting it 
to zero we obtain a set of equations that can be succinctly 
written in matrix form as 



[A + A(6)A b ] -C = A(a)A a • C. 



(12) 



C is a column vector with elements a's. Eq. 12 is a gener- 
alized eigenvalue equation which should yield exactly N 
real solutions. But the matrices A a and A b have rank 
one due to their peculiar structures [17], so Eq. 12 has a 
unique solution. We solve Eq. 12 using the computational 



software MATLAB to obtain the eigen-vectors C and C 
for two different values of the logarithmic derivatives A (b) 
and A (6). Note that the values of the logarithmic deriva- 
tives X(£) and A(^) at a given point x — i can be chosen 
arbitrarily [17]. In this paper, we have chosen £ = b for 
our calculations. 

Thus, in region II, we obtain two linearly indepen- 
dent solutions ip2{x) — ^2^=1 CiXi(x) and ip 2 (x) — 

~!2iLi CiXi(x)- Next we use the continuity of the wave- 
functions and their first derivatives at the boundaries 
x = a and x = b. These yield a set of four simultaneous 
equations 



(13) 



N 




E 


CiXi(a), 


JV 




E 


CiX' t (a), 


i=l 




N 




E 


CiXi(b), 


JV 




E 

i=l 


CiXt(b), 



(14) 



(15) 



(16) 



which we solve to evaluate the coefficients a\, b\, CI3, and 
b 3 . We repeat the same steps to obtain the other set 
of coefficients 0,1,61,03, and 63. These coefficients and 
q's completely determine the two independent wave- 
functions 4j(x) and 4>(x) for the entire region. We can 
now write the general solution for the entire region as a 
linear combination of these wave- functions, i.e. ^(x) — 
Bip(x) + Btp(x) where B and B are two arbitrary con- 
stants. 

Next, we obtain expressions for the reflection and 
transmission coefficients. In region I, using Eqs. 2-5, 
the general solution can be written as 

*i(;r) = B^ x {x) + B^x), 

= {B ai + Bax) [e tklX + pe - lklX ] . (17) 

p is the amplitude for reflection and is given by p — 
(Bbi + Bbi)/(Bai + Bdi). Similarly for region III, 



* 3 (a0 = B^ z {x)+B^{x), 

Ba 3 + Bd 3 ] e lk3X + \Bb. 



Bb ? 



3 e 



■ (18) 



We consider the physical situation where we have an inci- 
dent wave of unit amplitude (which implies i?oi +Bdi = 
1) coming from the left that partially penetrates the bar- 
rier and gets transmitted to region III. As there is no 
reflected wave in region III, we have another constraint 
Bb 3 + Bb 3 = on the coefficients B and B. By solving 
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these equations for B and B we obtain the reflection and 
transmission coefficients. They are given by 



R=\ P \- 



T= ^ 



ai& 3 - b z ai 



ai& 3 - 6 3 ai 



(19) 



(20) 



A good check of the numerical accuracy of our calcu- 
lation is to calculate the average energy E^ v in region II 
(a < x < b) with the input value of the energy E. E w 
can be evaluated using the expression 



E m = 



E 



JV p 
i,j = l c i c j r ij 



where 



Pij 



Rij - 

J a 



r b r n 2 



dxxi(x)xj(x). 



(21) 

(22) 
(23) 



A. Truncated linear step potential 

The truncated linear step potential barrier is a prac- 
tical example of a potential encountered in planar wave- 
guides [19]. It is defined by the potential function 



V(x) 



0, 

V (B-x), 
0, 



which is shown in Fig. 1. 

1.2 i , r- 



X 

> 



x < a 
a < x < b 
x > b 



V (B-x) 




(24) 



III. APPLICATIONS AND RESULTS 

In order to test the effectiveness of the variational 
i?-matrix method we have chosen five different one- 
dimensional potentials: truncated linear step, truncated 
exponential step, truncated parabolic, bell-shaped, and 
Eckart. We consider particles coming from region I with 
total energy E (E > Vi,Vs) striking the potential bar- 
rier. According to quantum mechanics, there is a finite 
probability for the particles to be transmitted to region 
III even if the particles have energy less than the max- 
imum height of the potential barrier. Particles that are 
not transmitted to region III are reflected back to region 
I. We use the variational method to calculate the reflec- 
tion and transmission coefficients as a function of particle 
energies for each potential barrier. 

For all cases we have checked the accuracy of our cal- 
culations by comparing the average energies calculated 
from Eq. 21 with the actual input energies. We also 
have compared our results with exact analytical results 
(ANA) and also with other numerical methods, such as 
numerical matrix (MAT) [9], transfer matrix (TM) [10], 
Wentzel-Kramos-Brillouin (WKB) [2,8], modified con- 
ventional WKB (MWKB) [7], modified Airy functions 
(MAF) [2], and improved MAF (MMAF) [18]. 



FIG. 1. Truncated linear step potential with Vb = 0.5 
hartree, B = 3 bohr, and boundaries at x = 1 bohr and 
x = 3 bohr. 

For the truncated linear step potential the Schrodingcr 
equation can be solved analytically to obtain the wave- 
function [20]. The exact solution for the inner region 
can be written in terms of the Airy functions Ai(x) and 
Bi(x) [21], 

4> 2 (x) = a 2 Ai (V /3 [B-x- e]) 

+ b 2 Bi (V/ 3 [B-x- e]) , (25) 

where a — 2mVo/h 2 1 e — E/Vq, and a 2 and b 2 are coef- 
ficients to be determined from the boundary conditions. 
One can then obtain the exact expressions for the reflec- 
tion (R e ) and tunneling (T e ) coefficients (see Appendix 
V). 

The reflection and tunneling coefficients calculated 
from our variational i?-matrix and analytical methods 
are plotted in Fig. 2 for different values of particle en- 
ergy E (in units of hartree) . As E increases it becomes 
easier for the particles to tunnel through the barrier - 
thus the tunneling coefficient increases (and the reflec- 
tion coefficient decreases as R + T = 1). For E w 1.1 
hartree, the reflection and tunneling probability curves 
cross each other. At this value of energy the probabili- 
ties of reflection and transmission are each equal to 50%. 
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When the energy of the particles are just over the maxi- 
mum height, VqB of the truncated linear step potential, 
i.e. E w 1.5 hartree, we would expect classically to see 
100% transmission. But from Fig. 2 we find the trans- 
mission probability to be about 77%. This reflection of 
particles for energies greater than the barrier height is 
also a quantum mechanical effect just as the transmis- 
sion of particles for energies less than the barrier height 
is a quantum mechanical effect. 




FIG. 2. Transmission T (filled circles) and reflection R 
(filled squares) coefficients through the truncated linear step 
potential barrier as a function of particle energy E. The solid 
(long-dashed) lines are for the exact values of T e (-Re) ob- 
tained from the analytical expression. For the wave-functions 
ip2(x) and ip2{x) we have used A(fe) = 1, \(b) — 4, and 
Ki — {0.1 : 0.1 : 6.0}(bohr) _1 as our input parameters. 

As a check of the accuracy of our calculation we have 
calculated the average energy E av using Eq. 21 and com- 
pared with the actual input energy E and have found 
that errors are in the range 0.03% - 1.75%. The higher 
errors occur for smaller values of the input energy, which 
suggests that more Ki's are required. 

In Table I we have compared the tunneling coefficients 
for different values of B obtained from our variational 
method (VRM) with results from other methods. The 
reason for choosing different values of B for the truncated 
linear step potential is to compare with existing data 
from other numerical methods [10,18]. It is known that 
the semi-classical WKB approximation fails close to the 
turning points (i.e. whenever E — V(x) ~ 0) - so those 
points are not taken into account in the WKB calcula- 
tion. Hence, changes in the B values have no effect on the 
transmission coefficient. In the modified WKB (MWKB) 
approximation the turning points are included and thus 
the transmission coefficient changes with changes in the 
B values. However the results are 5% - 14% higher than 



TABLE I. Comparison of the variation of tunneling coef- 
ficients for different values of B with different methods for 
truncated linear step potential. 



B 


VRM 


ANA 


TM 


MAT 


WKB 


MWKB 


MAF 


3 


0.6079 


0.6077 


0.6077 


0.6075 


0.4670 


0.6962 


0.6078 


6 


0.4819 


0.4818 


0.4818 


0.4815 


0.4670 


0.5163 


0.4856 


8 


0.4374 


0.4370 


0.4370 


0.4374 


0.4670 


0.4522 


0.4377 


15 


0.3424 


0.3423 


0.3423 


0.3416 


0.4670 


0.3353 


0.3431 


20 


0.3053 


0.3052 


0.3052 


0.3048 


0.4670 


0.2917 


0.3050 



the expected results. The variational results arc in ex- 
cellent agreement with the exact results and results from 
the numerical transfer matrix (TM), matrix (MAT), and 
modified Airy functions (MAF) methods. A good agree- 
ment of MAF results with the exact results is expected, 
as for the truncated linear step potential the Airy func- 
tions are the exact solutions of the Schrodingcr equation. 



B. Truncated exponential step potential 

The truncated exponential step function, shown in 
Fig. 3, is another example encountered in planar wave- 
guides [19]. The potential is characterized by the function 



V(x) 




x < a 
a < x < b 
x > b 



(26) 




FIG. 3. The truncated exponential step potential with 
Vo — 0.5 hartree. The potential is truncated at a = 1 and 
b = 8 bohr. 

For this potential an analytical solution of the 
Schrodinger equation exists, which can be written in 
terms of the modified Bessel functions of first and sec- 
ond kinds J(i/, z) and K(v, z) [21], 
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ip 2 {x) = a 2 K (2ik 1 2a 1/2 e- (x - a)/2 ^ 

+ b 2 I (-2ik, 2a 1 ' 2 e- {x - a ^ 2 ^ , (27) 



where k = y2mE/h and a — 2mVo/h . Using these 
wave- functions we have calculated the exact values of 
T e and R e (see Appendix V). Fig. 4 shows the results 
of our variational calculation and the analytical calcu- 
lation. The reflection and tunneling probability curves 
cross each other for E w 0.13 hartree. Over the max- 
imum height, Vo of the potential, i.e. E ss 0.5 hartree 
we find from Fig. 4, the transmission probability to be 
about 92%. Compared to the truncated linear step po- 
tential the transmission coefficient increases at a faster 
rate for the truncated exponential step function. This is 
expected as the slope of the truncated exponential poten- 
tial decreases as — (Vo/a) exp(— x/a) whereas the slope of 
the truncated linear step potential remains constant at 
-V Q . 




0.4 0.6 
Energy E 



FIG. 4. Transmission T (filled circles) and reflection R 
(filled squares) coefficients through the truncated exponen- 
tial step potential barrier. Solid (long-dashed) lines are for 
the exact values of T e (-Re) obtained from the analytical ex- 
pression. Different input parameters are A (b) = 2, A (6) = 8, 
and m = {0.1 : 0.1 : 6.0}. 



TABLE II. Comparison of the variation of tunneling coef- 
ficients for different values of E /Vo with different methods for 
truncated exponential step potential. 



E/V 


VRM 


ANA 


TM 


MAT 


WKB 


MWKB 


MAF 


0.25 


0.4812 


0.4789 


0.4789 


0.4788 


0.2247 


0.3892 


0.4865 


0.50 


0.7566 


0.7549 


0.7549 


0.7549 


0.4221 


0.8443 


0.7601 


0.75 


0.8712 


0.8702 


0.8702 


0.8695 


0.5693 


0.9861 


0.8733 



ational i?-matrix method are within 0.56% of the ana- 
lytical results (ANA); however, one should note that the 
agreement is improved as E increases. Other methods 
such as WKB and MWKB again fail to provide satis- 
factory results. Although MAF results are much more 
accurate than the WKB and MWKB results but it is 
less accurate (differs from exact results by 0.35% - 1.5%) 
than our variational results. 



C. Truncated parabolic potential 

The truncated parabolic potential illustrated in Fig. 5 
is of considerable interest in fiber and integrated optical 
waveguides [19]. The profile of the potential is given by, 

{0, x < a 

Vo [B 2 -{x- x ) 2 ] , a < x < b (28) 
0. x>b 



0.6 



0.4 



V(x)=V [B 2 -(x-x ) 2 ] 



x 
> 



0.2 




We have compared the calculated E av values with the 
input E values to check the accuracy of our calculation 
and have found that the error is between 0.013% - 2.1% of 
the input energy values. This again shows that a denser 
grid of Ki's is needed for smaller E values. 

Our variational results (VRM) are compared to other 
methods in Table II [10,18]. Here we find that the nu- 
merical transfer matrix (TM) and matrix (MAT) meth- 
ods give the same results as the analytical results. An 
inspection of Table II shows that the results of our vari- 



FIG. 5. Parabolic potential truncated at x — xo ± B bohr. 
We have chosen Vo = 0.5 hartree, xo — 2 bohr and B — 1 
bohr for our calculation. 

In Fig. 6 we have plotted our variational i?-matrix re- 
sults for the transmission and reflection coefficients. Here 
the particle with E 0.27 hartree have equal reflection 
and transmission probabilities. The transmission prob- 
ability is 74% when the particle energy just exceeds the 
maximum height VoB 2 of the barrier. 
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Energy E 

FIG. 6. Transmission T (filled circles) and reflection R 
(filled squares) coefficients through the truncated parabolic 
potential barrier. The exact values for T e and R e are not 
known in this case. The lines are drawn only to guide the 
eyes. Input parameters used are X(b) — 2, X(b) — 3, and 
m = {0.1,0.5,0.9,1.3,2.0,2.4,3.0,3.4,4.0,4.4,5.0,5.4,6.0} 
(bohr)- 1 . 

TABLE III. Comparison of the variation of tunneling co- 
efficients for different values of E/Vo with different methods 
for truncated parabolic potential. 



E/V 


VRM 


TM 


MAT 


WKB 


MAF 


MMAF 


0.10 


0.1122 


0.1124 


0.1122 


0.0575 


0.0390 


0.0971 


0.20 


0.2139 


0.2141 


0.2141 


0.0778 


0.0506 


0.1859 


0.50 


0.4599 


0.4604 


0.4603 


0.1878 


0.1003 


0.3981 



When compared to our E av values with the input E 
values we obtained a satisfactory agreement as the errors 
are within 2%. Our variational results are compared with 
results from other numerical methods for three different 
values of E/V in Table III [10,18]. The exact analyt- 
ical results for T e and R e are unavailable in this case. 
This comparison show that our results from variational 
i?-matrix method is in excellent agreement with trans- 
fer matrix and matrix methods. In this case also, WKB, 
MAF, and MMAF does not provide any satisfactory re- 
sults. 



where x is a constant. 
0.6 




FIG. 7. Bell-shaped potential with Vo = 2 hartree and 
= 5 bohr. Inner region is defined between x = 1 bohr 
and x = 9 bohr. Potentials V\ in region I and Vs in region III 
are assumed to be zero. 

An exact solution of the wave-function of the particle 
in region II for this potential barrier is known and can 
be written in terms of the hypergeometric functions [20] . 
Transmission and reflection coefficients can then be cal- 
culated. The exact tunneling coefficient is given by (for 
8Vq > 1, which is our case) the closed form expression, 



T e = 



sinh 2 (7rA:) 



sinli (Trfc) +coslr (tt/3/2) 



(30) 



where k = \ 2mE/ti and (5 



8mV /h 2 - 1. In this 



case our calculated E av values differ from the input E 
values by 0.0065% - 0.694%. For the truncated bell- 
shaped potential we are unable to compare with other 
numerical methods as we have not found any existing 
numerical results. Fig. 8 illustrates our variational results 
and the exact results. At larger values of the energy E 
our values of T and R deviate slightly from the exact 
values which again suggests that a better and optimized 
choice of n{' is necessary to obtain better accuracy. From 
Fig. 8 we find that the transmission coefficient is about 
62% when E is over the peak of the barrier. 



E. Eckart's potential 



D. Truncated bell-shaped potential 

Our next example is the smooth bell shaped potential 
defined by 



V(x) 



Vo 



cosh (x — xo) 



(29) 



Eckart introduced a smoothly varying potential func- 
tion defined as 



V( X ) = - 



1 + e x 



+ B- 



(1 + e x 



(31) 



where A and B are two dimensionless arbitrary parame- 
ters [22]. The shape of the potential is shown in Fig. 9. 
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Energy E 

FIG. 8. Transmission T (filled circles) and reflection R 
(filled squares) coefficients through the bell-shaped potential 
barrier. Solid and long-dashed lines are for the exact values 
of T e and R e . Input parameters used are A(fo) = 1, A(b) = 9, 
and K t = {0.1 : 0.2 : 3.0}. 

The Schrodingcr equation for this potential leads to a 
differential equation of hypergeometric type, whose solu- 
tions can be written down in terms of the hypergeometric 
functions from which the exact coefficients of reflection 
and transmission can be extracted [1,22]. The exact ex- 
pression of R e is given by (for B > 1/4) 

_ cosh(2^(fc-/?))+cosh(27n5) 
e cosh(27r(fc + $)) + cosh(27rS) ' 1 ' 



where k = yj 2mE/h 2 ,[3 = Vk 2 - A, and S = ^JB - 1/4. 
The coefficient T e is obtained using the relation T e = 
1 — R e . The computed values of T, R and the exact val- 
ues T e ,R e are shown in Fig. 10. Our variational results 
for T and R are in good agreement with the exact re- 
sults. From Fig. 10 we find that almost all the particles 
are transmitted to region III when their energies are just 
above the maximum height of the barrier, which is at 2.23 
hartree. For large values of E there are some discrepan- 
cies between the two results which could be improved 
by enlarging and optimizing the set of m. Note that in 
this case also we have not found any data using other 
numerical techniques for comparison. For this case, the 
errors for i? av compared to the input E values are within 
0.078% - 1.43%. 



IV. CONCLUSION 

We have applied the variational i?-matrix method to 
calculate the reflection and tunneling probabilities of par- 



2.4 




FIG. 9. Eckart potential with A = 1, B = 8, and x = 8 
bohr. Inner region is defined between x = 2 bohr and x = 13 
bohr. Potentials Vi and V3 are approximately constants 
and their values are evaluated at the boundary values, i.e. 
V\ — V{x)\ x= \ and V3 = V{x)\ x= \z. The particle energies E 
are taken to be larger than Vi and V3. 

tides tunneling through five one-dimensional potential 
profiles - truncated linear step, truncated exponential 
step, truncated parabolic, bell shaped, and Eckart. We 
have compared our results with results from analyti- 
cal and other numerical methods and found satisfactory 
agreement with the analytical results and some of the 
numerical results obtained using numerical matrix and 
transfer matrix techniques. Our results can be improved 
by a better choice of a larger and optimized basis set. We 
have checked this by enlarging the basis sets for the trun- 
cated linear and exponential potential steps. Although 
semi-classical WKB method is a simple method but it 
is restricted to slowly varying potentials that are con- 
tinuous. For potentials with discontinuities WKB ap- 
proximation fails. In modified WKB approximation this 
deficiency is taken into consideration, but still it fails 
to provide satisfactory results for any general potential. 
Modified Airy function and improved MAF methods are 
suitable only for certain potentials such as the truncated 
linear step, truncated parabolic, and truncated quartic 
potentials. Although the matrix method and the transfer 
matrix technique yield satisfactory results but they are 
complicated. On the other hand the variational i?-matrix 
method is simple, general, easy to implement, and applies 
to any potential shape. Also it is non-iterative and thus 
computationally efficient. This method can be applied in 
scattering theory and also in quantum tunneling in two 
and three dimensions. 

The authors thank S. T. Powell and A. Lahamer for 
helpful discussions and suggestions and gratefully ac- 
knowledge the support by the Undergraduate Research 
and Creative Projects Program of Berea College. 
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FIG. 10. Transmission T (filled circles) and reflection R 
(filled squares) coefficients for the Eckart potential. Solid and 
long-dashed lines are for the exact values of T e and R e ob- 
tained from the analytical expression in Eq. 32. Input param- 
eters used for this calculation are A(fe) = 2, X(b) = 13, and 
Ki = {0.001 : 0.2 : 3.0}. 



[1] M. Razavy, Quantum Theory of Tunneling, World Scien- 
tific, NJ, 2003. 

[2] A. K. Ghatak and S. Lokanathan, Quantum Mechanics: 
Theory and Applications, 3rd. Ed., New Delhi; McMillan, 
1984. 

[3] A. Casher and F. Englert, Class. Quantum Grav. 10 
2479, (1993). 

[4] D. S. Chemla and A. Pinczuk, Eds. IEEE J. Quan- 
tum Electron., Special Issue on Semiconductor Quantum 
Wells and Superlattices: Physics and Applications, QE- 
22, 1986. 

[5] S. Datta, Electronic Transport in Mesoscopic Systems, 

Cambridge Univ. Press, N. Y, 1995. 
[6] D. D. Awschalom, D. Loss, and N. Samarth (Eds.), 

Semiconductor Spintronics and Quantum Computation, 

Springer, N. Y., 2002. 
[7] J. D. Love and C. Winkler, J. Opt. Soc. Am. 67, 1627 

(1977). 

[8] A. K. Ghatak, R. L. Gallawa, and I. C. Gayal, IEEE J. 

Quantum Electron. 28, 400 (1992). 
[9] A. K. Ghatak, K. Thyagarajan, and M. R. Shenoy, 

J. Lightwave Technol. 5, 660 (1987); A. K. Ghatak, 

K. Thyagarajan, and M. R. Shenoy, IEEE J. Quantum 

Electron. 24, 1524 (1988). 
[10] A. Zhang, Z. Cao, Q. Shen, X. Dou, and Y. Chen, J. 

Phys. A: Math. Gen. 33, 5449 (2000). 



[11] A. M. Lane and R. G. Thomas, Rev. Mod. Phys. 30, 
(1958). 

[12] E. P. Wigner and L. Eisenbud, Phys. Rev. 72, 29 (1947). 

[13] L. Smrcka, Microstruct. 8, 221 (1990). 

[14] W. Kohn, Phys. Rev. 74, 1763 (1948). 

[15] R. K. Nesbet, Variational Methods in Electron-Atom 

Scattering Theory, Plenum, N. Y, 1980. 
[16] H. L. Rouzo and G. Raseev, Phys. Rev. A 29, 1214 

(1984). 

[17] H. L. Rouzo, Am. J. Phys. 71, 273 (2003). 

[18] S. Roy, A. K. Ghatak, I. C. Goyal, and R. L. Gallawa, 
IEEE J. Quantum Electron. 29, 340 (1993). 

[19] M. S. Sodha and A. K. Ghatak Inhomogeneous Optical 
Waveguides, Plenum, New York, 1977. 

[20] L. D. Landau and E .M. Lifshitz, Quantum Mechanics: 
Non-relativistic Theory, Addison- Wesley, MA., 1958. 

[21] M. Abramowitz and I. A. Stegun Handbook of Mathemat- 
ical Functions, Dover Publications, New York, 1970. 

[22] C. Eckart, Phys. Rev. 35, 1303 (1930). 



V. APPENDIX 

We assume that for a potential V(x) the Schrodinger 
equation (Eq. 1) can be solved exactly for the inner region 
defined by a < x < b and the solution is given in terms 
of two functions f(x) and g(x) as 



ip2{x) = a 2 f{x) + b 2 g{x). 



(33) 



For example, for the truncated linear step potential the 
functions f(x) and g(x) are the Airy functions Ai(x) 
and Bi{x) and for the truncated exponential step poten- 
tial they are the modified Bessel functions K (y, z) and 
I{y,z). We now match the wave functions tpi(x),ip2(x) 
(■tp2(x),tp 3 (x)) and their derivatives at x — a (at x = b) 
and solve the set of simultaneous equations to obtain the 
coefficients dj and b,'s (i = 1,2,3). Using these coef- 
ficients an exact expression for the tunneling coefficient 
can be obtained which is given by, 

(34) 

where the numerator N and the denominator D are, 
N=[f(b)g'(b)-g(b)f'(b)}, 
D = 





N 




~D 



g'(b) - ik 3 g(b) 
f'(b) - ihf(b) 



1 - jr/<°\ 
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